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ABSTRACT 

We have obtained 18 new high-resolution spectropolarimetric observations of the 
B2Vp star a Ori E with both the Narval and ESPaDOnS spectropolarimeters. The 
aim of these observations is to test, with modern data, the assumptions of the Rigidly 
Rotating Magnetosphere (RRM) model of Townsend & Owocki (2005), applied to the 
specific case of a Ori E by Townsend et al. (2005). This model includes a substantially 
offset dipole magnetic field configuration, and approximately reproduces previous ob- 
servational variations in longitudinal field strength, photometric brightness, and Ha 
emission. We analyze new spectroscopy, including H I, He I, C II, Si in and Fe III 
lines, confirming the diversity of variability in photospheric lines, as well as the dou- 
ble S-wave variation of circumstellar hydrogen. Using the multiline analysis method 
of Least-Squares Deconvolution (LSD), new, more precise longitudinal magnetic field 
measurements reveal a substantial variance between the shapes of the observed and 
RRM model time- varying field. The phase resolved Stokes V profiles of Hei 5876 A and 
6678 A lines are fit poorly by synthetic profiles computed from the magnetic topol- 
ogy assumed by Townsend et al. (2005). These results challenge the offset dipole field 
configuration assumed in the application of the RRM model to a Ori E, and indi- 
cate that future models of its magnetic field should also include complex, higher-order 
components. 

Key words: stars: magnetic fields - stars: rotation - stars: early-type - stars: circum- 
stellar matter - stars: individual: HD 37479 - techniques: spectroscopic 



* Based on observations obtained using the Narval spectropo- 
larimeter at the Observatoire du Pic du Midi (France), which 
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of the Centre National de la Recherche Scientifique of France, and 
the University of Hawaii, and observations collected at the Eu- 
ropean Organisation for Astronomical Research in the Southern 
Hemisphere, Chile, under program ID 62.H-0319 and during the 
FEROS Commissioning run II. 
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1 INTRODUCTION 

Hot, massive stars are not expected to to harbor magnetic 
fields due to the absence of strong convection zones in their 
outer envelopes. However, peculiar A and B (Ap/Bp) stars 
are known to possess strong, organized magnetic fields with 
strengths up to tens of kG. The hottest of these stars, the 
He-strong stars, are main sequence B2 stars with enhanced 
and often variable helium abundance. 

The helium-strong B2Vp star a Orionis E (HD 37479) 
is one of the most well-studied variable massive stars. This 
bright star is characterized by rapid rotation, a strong mag- 
netic field, and variable circumstellar Ha emission. Obser- 
vations show modulation according to the 1.19 day rotation 
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period in longitudinal magnetic field (Landstreet & Borra 
1978), Ha emission (Walborn 1974), He line strength (Ped- 
ersen & Thomsen 1974), photometry (Hesser et al. 1976), 
UV line strength (Smith & Groote 2001), 6 cm radio emis- 
sion (Leone & Umana 1993), and linear polarization (Kemp 
& Herman 1977). a Ori E has long been known as the pro- 
totypical magnetic Bp star. 

Over three decades ago, Landstreet & Borra (1978) dis- 
covered the strong magnetic field of a Ori E, the first detec- 
tion in a helium-strong star. Using the circular polarization 
in the wings of the H/3 line, the eight longitudinal magnetic 
field (B^) measurements, when phased with the photometric 
ephemeris of Hesser et al. (1976), form a curve that appears 
nearly sinusoidal, varying from -2.2 kG to 2.8 kG. Assum- 
ing a centered dipole magnetic field, the authors estimated 
a stellar surface field strength at the magnetic poles of ~10 
kG. As a follow-up study to this seminal paper, Bohlender et 
al. (1987) observed ten helium-strong stars, presenting eight 
H/3 and four He I 5876 A circular polarization measurements 
of a Ori E. The additional data confirms the periodicity and 
field amplitude of Landstreet & Borra (1978). 

Walborn & Hesser (1976) suggested that this star was 
in a mass-transfer binary system with an unseen compan- 
ion. However, the authors also discussed a Ori E as a possi- 
ble application of the single-star oblique rotator model pro- 
posed by Stibbs (1950). With the landmark discovery of the 
star's field, Landstreet & Borra (1978) proposed a frame- 
work which depicts a Ori E as a magnetic oblique rotator 
with plasma trapped in a magnetosphere. Theoretical efforts 
to establish a viable model to describe the physical phe- 
nomena included work by Nakajima (1985), Shore & Brown 
(1990), Bolton (1994), Short & Bolton (1994), Preuss et al. 
(2004), and most recently, Townsend & Owocki (2005). 

The Rigidly Rotating Magnetosphere (RRM) model 
was developed by Townsend & Owocki (2005) to analytically 
describe the circumstellar plasma structure of a rapidly ro- 
tating star in which the magnetic field overpowers the stellar 
wind, causing plasma to become trapped in a magnetosphere 
that co-rotates with the star. The strong magnetic field spins 
up and confines the wind plasma, keeping it rigidly rotat- 
ing well beyond the Kepler co-rotation radius, but also held 
down against the net outward gravito-centrifugal force. The 
geometry of the magnetosphere depends on the obliquity 
angle j3 and the angular velocity SI of the star (Townsend 
2008). The observed variability also depends on the star's 
inclination angle, i. 

Townsend et al. (2005) applied the RRM model to the 
specific case of a Ori E, adopting the parameters i = 75°, 
O = 0.5S1 C , and /3 — 55°. Q, c is the velocity at which the 
surface centrifugal force at the equator balances with grav- 
ity. The benefit of choosing such a well-studied star is the 
ability to robustly confront the model with a variety of ob- 
servations. Model parameter determinations were strongly 
driven by observational constraints. Motivated by the un- 
equal primary and secondary minima in the observed pho- 
tometric light curve, as well as the unequal strength of the 
blue and red-shifted Ha emission, the authors defined an 
offset vector, a, to describe the displacement of the dipole 
magnetic center from the stellar center. In units of the polar 
radius R p , the dipole magnetic field was significantly offset 
to a = (-0.041, 0.30, 0.029) (corresponding to a 0.3 R p dis- 
placement perpendicular to the rotating and magnetic axes 



and a —0.5 R p displacement along the magnetic axis). In 
confirmation of original suggestions by Landstreet & Borra 
(1978), the RRM model adopted for a Ori E produces two 
regions of enhanced plasma density (clouds), located at the 
intersections of the magnetic and rotational equators. A sub- 
stantial success of the model was its ability to qualitatively 
reproduce the observed Ha variations of the FEROS spectra 
previously analyzed by Reiners et al. (2000). The model pho- 
tometric light curve agrees with the observed eclipse depths 
and relative timing of the Hesser et al. (1977) Stromgren u- 
band data. The measurements of Landstreet & Borra (1978) 
and Bohlender et al. (1987) correspond reasonably well to 
the computed model longitudinal field curve. However, dif- 
ferences do exist between the model and observed data, such 
as the quantitative properties of the Ha variations. The 
excess photometric brightness at rotational phase ~0.6 is 
likely a result of photospheric features not considered by 
this purely circumstellar model, as shown by Krtifika et al. 
(2007) for HD 37776. 

The goal of this project is a re-evaluation of the RRM 
model using modern data to test the original assumptions 
of Townsend et al. (2005), and to resolve the remaining dis- 
crepancies between the data and the model. In this first 
paper, we present new spectropolarimetric observations of 
a Ori E. Section 2 describes the observations and data re- 
duction. Section 3 explains the magnetometry analysis and 
results. Section 4 explores the spectroscopic properties. We 
discuss in Section 5 the RRM model's ability to accurately 
describe a Ori E in light of these new observations. Section 
6 presents a summary of this paper, as well as future work. 



2 OBSERVATIONS 

We obtained a total of 18 high-resolution (R=65000) broad- 
band (370-1040 nm) spectra of a Ori E. Sixteen spectra 
were obtained in November 2007 with the Narval spectropo- 
larimeter on the 2.2m Bernard Lyot telescope (TBL) at 
the Pic du Midi Observatory in France. The remaining two 
spectra were obtained in February 2009 with the spectropo- 
larimeter ESPaDOnS on the 3.6-m Canada-France-Hawaii 
Telescope (CFHT), as part of the Magnetism in Massive 
Stars (MiMeS) Large Program (Wade et al. 2010). The ob- 
servation log is reported in Table [T] Reduction was per- 
formed at both observatories with the Libre-ESpRIT pack- 
age (Donati et al. 1997), which yields the Stokes I (inten- 
sity) and Stokes V (circular polarization) spectrum, as well 
as the null spectrum (N), which diagnoses any spurious con- 
tribution to the polarization measurement. Four consecutive 
sub-exposures are combined using double ratios to produce 
one polarization spectrum. Each resultant Stokes I spectrum 
was normalized by individual order using a polynomial fit to 
continuum points. Each Stokes V spectrum was normalized 
by the same polynomial continuum fit, producing a V/I c 
spectrum. 

Additionally, 29 FEROS echelle spectra were used in 
the spectroscopic analysis of a Ori E (Kaufer et al. 1999; 
See Section 4). These data were taken in November and 
December 1998, partially during the commissioning of the 
FEROS spectrograph on the ESO 1.52-m telescope at La 
Silla. These observations are further described by Reiners et 
al. (2000). 
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Table 1. Log of spectropolarimetric observations and magnetic field measurements. All but the last two rows of data were taken with 
the Narval spectrpolarimeter. The remaining two observations were taken with the ESPaDOnS spcctropolarimeter. Column 1 gives the 
heliocentric Julian date of mid-observation. Column 2 the rotational phase according to the adopted ephcmcris of Townsend et al. (2010). 
Column 3 the total exposure time of the observation. Column 4 the peak signal to noise ratio per 1.8 km s — 1 per spectral pixel in the 
reduced spectrum. Column 5 the signal to noise ratio in the LSD profile. Columns 6, 7 and 8 the longitudinal magnetic field computed 
from individual line Stokes V profiles for H/3, He I 6678 A, and He I 5876 A, along with their associated formal errors. Columns 9 and 
10 the longitudinal magnetic field computed from LSD mean Stokes V and N profiles, along with their associated formal errors. 
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3 MAGNETIC FIELD MEASUREMENTS 

3.1 Least Squares Deconvolution 

We used the multiline analysis method of Least-Squares De- 
convolution (LSD; Donati et al. 1997) to produce mean 
Stokes / and V profiles from our spectra for calculation 
of high-precision longitudinal magnetic field measurements. 
Silvester et al. (2009) describe the details of the LSD method 
as applied to B stars. To determine appropriate parameters 
for the LSD mask, we studied the effect of various combi- 
nations of T e ff, loggr, line depth cutoff level, spectral ranges, 
and chemical elements on the LSD results. 

Starting with Vienna Atomic Line Database (VALD; 
Kupka et al. 1999) "extract stellar" line lists, we produced 
synthetic LTE spectra using the spectral synthesis code 
Synth3 described by Kochukhov (2007). For consistency, 
each synthetic spectrum was compared to the same observa- 
tion, taken at phase 0.672 during minimum helium absorp- 
tion. The spectra used initial parameters of T c g — 23000 K, 
loggr = 4.0, as suggested by Cidale et al. (2007), although 
spectra were also computed by varying both temperature 
(by ± 3000K) and log g (by ± 0.5). The abundance parame- 
ters of the synthetic spectra were adjusted until good visual 
agreement was reached with the observed spectra. Although 
the He lines in the red part of the spectrum are affected by 
Stark broadening and NLTE effects, the bluer helium lines 
were best fit by a helium abundance of log(NHo/Nrot) = 
—0.1. Although lines of certain metals may be variable, so- 
lar abundances fit well for a majority of rotational phases. 
At this stage, we set the chemical abundances and proceeded 
to determination of an LSD line mask. 

We next created line masks consisting of spectral line 
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Figure 1. LSD Stokes I (bottom), null N (middle), and Stokes V 
(top) profiles of o Ori E at rotational phase 0.672. The profiles 
are expanded by the indicated factor. The error bars on the null 
and Stokes V data are propagated from the formal error bars 
associated with each pixel. For visualization purposes, the null 
and Stokes V have been shifted upwards from 0.0. A clear Zeeman 
signature is detected in the Stokes V profile, while the null profile 
shows no signal. 



wavelengths, line depths, and atomic data for each transi- 
tion contributing to the spectrum. The normalized obser- 
vational spectra were then processed using this line mask 
in the LSD procedure. The output of LSD provides mean 
Stokes I and V profiles, mean null profiles, as well as statis- 
tics related to the procedure fit and the output profiles. The 
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LSD procedure propagates the formal error associated with 
each pixel throughout the deconvolution process. An exam- 
ple of the LSD profiles for the final line mask is displayed 
for phase 0.672 in Fig. [T] We adopted a velocity bin of 2.6 
km s _1 for the calculation of profiles. The velocity range 
was set from —243 km s _1 to +343 km s _1 to include the 
line itself and check for any signal that may be outside the 
profile. 

To test the sensitivity of the LSD line mask, we con- 
ducted an experiment to explore the impact of varying cer- 
tain parameters. A large set of masks were created and run 
through the LSD procedure to view any differences. Inclu- 
sion or exclusion of certain lines or elements is an important 
effect on the quality of the LSD profiles. When we excluded 
helium lines from our line mask, the resulting profiles were 
extremely noisy, as a result of the weakness of metal lines. In- 
cluding lines located in the spectral regions overlapping the 
higher- order B aimer series (A < 4000 A) or regions domi- 
nated by telluric lines (A > 8000 A) also affected the results 
negatively. Varying T e ff by ± 3000 K and log<? by ± 0.5 
had little effect on the profiles. A line depth limit of greater 
than 10% of the continuum produced the cleanest LSD pro- 
files, however changing the value did not significantly worsen 
the results. At the conclusion of the investigation, the best 
overall line mask was chosen for the spectrum of a Ori E. 

The final mask we selected begins with a line list for 
T e g = 23000 K, log g — 4.0, and solar abundances except 
for helium. The helium abundance was set to log(NH e /Nxot) 
= —0.1. A total of 332 lines contributed to the line mask, 
and most helium and metal lines in the range 405-800 nm 
were included. Lines blended with hydrogen Balmer lines 
within the included range were removed, as hydrogen lines 
are not included in the mask. The wavelength range was 
chosen in an attempt to exclude regions with many hydro- 
gen lines in the blue part of the spectrum and strong telluric 
contamination in the red part of the spectrum. To be clear, 
this study is not meant as a precise stellar parameter de- 
termination. We have not modeled in detail the spectrum 
and are not attempting to make any assertions regarding 
effective temperature or gravity of this star. The process de- 
tailed here simply attempts to determine the overall best 
line mask to extract clean Stokes I and V profiles from the 
LSD procedure. 

Lack of signal in the null profiles indicates that there are 
no important spurious contributions to the Stokes V profiles. 
Each separate LSD Stokes V profile produces a definite de- 
tection (detection probability > 99.999%) according to the 
criteria described by Donati et al. (1997). Fig. [^presents the 
Stokes / and V profiles for both the He I 6678 A line, as well 
as the LSD mean I and V profiles, for each rotational phase. 
The helium and LSD Stokes V profiles agree quite well, as 
expected as helium lines are the strongest lines contributing 
to the LSD mask. 



3.2 Longitudinal Magnetic Field 

Historically, the longitudinal magnetic field measurements 
(Bf) of a Ori E were computed from the circular polarization 
signatures in the wings of the H/3 and He I 5876 A lines (e.g. 
Landstreet 1980, 1982; Borra & Landstreet 1977). We utilize 
the circular polarization and intensity spectra of H/3, He I 
6678 A, and He I 5876 A to determine the B< value at each 



phase. These values are reported in Table [T] along with la 
uncertainties. We calculate the longitudinal magnetic field 
from the first moment of Stokes V: 



B e = -2.14 x 10 1 



J vV(v)dv 



Xgc J[I C - I(v)]dv 



G. 



(1) 



where g is the Lande factor and A is the wavelength of the 
specified line (Donati et al. 1997; Wade et al. 2000). The in- 
tegration range (in velocity, v relative to the line center-of- 
gravity) employed in calculating the longitudinal magnetic 
field was from —215 km s _1 to +215 km s _1 for helium lines 
and from -500 km s _1 to +500 km s" 1 for H/3. These ve- 
locity ranges were chosen to be sure to include the entire 
Stokes V profile, but not too wide as to artificially increase 
the error estimates by including excess continuum. The new 
measurements are plotted in Fig. [3] along with the corre- 
sponding historical data from Landstreet & Borra (1978) 
and Bohlender et al. (1987). The data are phased according 
to the ephemeris found by Townsend et al. (2010): 

JD = 2442778 d 819 + i d 1908229£ + l d 44 x 10~ 9 £ 2 . (2) 

This second-order ephemeris accounts for the increasing pe- 
riod of a Ori E due to rotational braking by the magnetic 
field. 

The historical H/3 measurements and new measure- 
ments from H/3 Stokes V profiles agree relatively well. There 
is better consistency with the He I lines. The new mea- 
surements have much smaller error bars by a factor of ~3, 
and the amplitude of the variation appears smaller than the 
historical data. The amplitude difference could be a con- 
sequence of the fact that the historical measurements were 
taken from the wings of the line, while our measurements 
are mainly in the core. A second order least-squares sine 
curve fit to the new H/3 measurements gives a reduced \ 2 
value of 1.54. The new helium data from He I 5876 A and 
6678A lines fit with a second order least-squares sine curve 
have a reduced x 2 °f 2.62. Fig [3] shows both the data and 
the corresponding fits. When examining measurements from 
two different elements, it is important to remember that lon- 
gitudinal magnetic field curves are affected not only by the 
line of sight component of the field, but also by the surface 
distribution of that element (Deutsch 1958; Pyper 1969). 
Bohlender et al. (1987) found little difference between the 
hydrogen and helium measurements as far as the shape of 
the curve and the magnetic extrema are concerned. We find 
a similar shape, but the H/3 measurements best fit has ex- 
trema ~ 150-250 G larger than the best fit to the helium 
data. 

From the mean LSD profiles, we used Eq. {TJ to calcu- 
late the longitudinal magnetic field, setting g and A to the 
S/N- weighted mean values for all the included lines. The 
integration range employed in calculating the longitudinal 
magnetic field was from —200 km s _1 to +200 km s _1 , with 
the same considerations as for the individual line measure- 
ments. The derived values, along with their la uncertainties, 
are reported in Table [1] The longitudinal magnetic field is 
observed to change from about —1.9 kG to +2.4 kG. Uncer- 
tainties for the longitudinal field values vary from 54 to 164 
G, decreasing the error bars by a factor of four compared 
to the historical data. The longitudinal magnetic field curve 
is shown in the bottom panel of Fig. [3] including the his- 
torical magnetic data. The agreement between new and his- 
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Figure 2. Stokes / and Stokes V profiles of a Ori E for the He I 6678 A line (left) and for the LSD fits (right). The profiles are expanded 
by the indicated factor and each separate observation is shifted upwards to its proper rotational phase calculated using the ephemeris 



in Eq. (f2j) . Clear Zeeman signatures are detected in each Stokes V profile. The He I 6678 
velocity, while the calculation of LSD profiles adopted a velocity bin of 2.6 km s _1 . 



profiles have been binned by 2 km s 



toric data, phased according to the Townsend et al. (2010) 
ephemeris, indicates long-term stability of the large-scale 
magnetic field. The photometric light curve structure shows 
similar stability, further indicating that the entire system, 
while not static, is not undergoing large-scale changes in 
magnetic and magnetospheric structure. Magnetic extrema 
occur at phases 0.2 and 0.8, while null longitudinal fields 



occur at phases 0.0 and 0.4. The shape of the magnetic field 
curve is quite clearly non-sinusoidal. A third order least- 
squares sine curve fits the data well (reduced x 2 — 2.26). 
The asymmetry of the longitudinal magnetic field variation 
may indicate a magnetic field that is not a simple oblique 
dipole. 
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Figure 3. Longitudinal magnetic field measurements for <r Ori 
E. The solid curve is the second order (H/3, He i) or third or- 
der (LSD) least-squares sine curve fit to the CFHT and Narval 
spectropolarimetric data. The data are phased according to the 
ephemeris of Townsend et al. (2010). Top: H/3 Stokes V mea- 
surements (filled circles), as well as measurements reported in 
Landstreet & Borra (1978) (squares) and Bohlender et al. (1987) 
(triangles) for H/3. Middle: He I 6678 A (filled circles) and He I 
5876 A (asterisks) Stokes V measurements, and helium line mea- 
surements from Bohlender et al. (1987) (diamonds). Bottom: 
LSD profile measurements (filled circles), as well as the data re- 
ported in Landstreet & Borra (1978) and Bohlender et al. (1987) 
(squares) 



Using the Lomb-Scargle periodogram (Lomb 1976; 
Scargle 1982), the historical and LSD derived magnetic data 
gives a period of 1.190842 ± 0.000004 days, where the re- 
ported error is the 3<r uncertainty. Our derived period is 
longer than the period of 1.19081 ± 0.00001 days Hesser et 
al. (1977) obtained from photometric data, but consistent 
with a period of 1.19084 ± 0.00001 days derived from He I 
4471 A equivalent width variations in Reiners et al. (2000). 
Our data span a baseline long enough that the period can be 
determined down to a third of a second. However, due to the 
magnetic braking reported by Townsend et al. (2010), our 
period determination from the longitudinal magnetic field 
data is strictly an average period over the time frame of 
the observations. In fact, when the current epoch of mag- 
netic data are phased with the mean derived period and the 
HJDo of Hesser et al (1976), the data do not define a null 
longitudinal magnetic field at phase 0.0, suggested by the 
RRM model. If the data are phased with the ephemeris of 
Townsend et al. (2010), as given in Eq. ((2]), the longitudinal 
field curve better reflects the expected phases of null fields 
and extrema. 



3.3 Summary 

We have derived longitudinal magnetic field measurements 
from Stokes V profiles of the H/3, Hel 5876 A, and He I 
6678 A lines. These new measurements agree relatively well 
with historical data, although the amplitude of the variation 
is slightly smaller for the new data. We also used the LSD 
method to produce clean, accurate mean Stokes /, Stokes V, 
and null profiles, as well as longitudinal magnetic field mea- 
surements for each of the 18 new spectropolarimetric obser- 
vations of a Ori E. The new measurements, when phased ac- 
cording to the ephemeris of Townsend et al. (2010), exhibit 
a non-sinusoidal periodicity (Fig. [3}. The overall shape of 
the curve suggests a deviation from a simple oblique dipole. 
Comparison with historical data indicates a long-term sta- 
bility of the magnetic field structure. In the discussion, we 
will compare the currently available observations with the 
offset dipole magnetic field configuration that Townsend et 
al. (2005) applied to the RRM model for a Ori E. 



4 SPECTROSCOPY 

The spectroscopic properties of a Ori E have been a popu- 
lar research focus for decades, since the discovery of strong 
helium lines by Berger (1956). With this new set of high- 
resolution spectra, our future goal is to model the pho- 
tospheric spectral variations in different elements using 
Doppler imaging, as well as magnetospheric variability. Here 
we examine the profiles and variability of the spectral lines 
of a Ori E in order to evaluate their qualitative properties, 
and identify lines appropriate for magnetic Doppler imag- 
ing (MDI; Piskunov & Kochukhov 2002). Due to the hot 
temperature of B2Vp stars, strong spectral lines in the op- 
tical spectrum are primarily H I and He I. We searched for 
weaker lines with high enough S/N and which showed vari- 
ability. Investigation of various lines of C II, N II, O II, Al ill, 
Si in, S II, and Fe ill support the idea that the variability of 
a Ori E extends to both circumstellar and photospheric en- 
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Figure 5. Dynamic spectral plots for Ho and H/3, phased using the ephemeris of Townsend et al. (2010). Spectra are plotted minus a 
synthetic photosphcric spectrum. Dashed vertical lines mark the stellar rotational velocity. Dashed curves represent the motion of the 
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vironments. Each of the lines exhibits variability, although 
some on a more pronounced scale than others. 

In a study of optical line variability, Reiners et al. (2000) 
used the FEROS commissioning spectra of a Ori E (Kaufer 
et al. 1999) to study and model photospheric lines of He I, 
C II, and Si III. The authors apply a model with a dipole 
magnetic field and two circular regions of variable abun- 
dance and size at the magnetic poles. Reiners et al. (2000) 
conclude that for a Ori E the magnetic polar regions (each 
with a radius of 60°) are overabundant in helium (log(He/H) 
= 0.0) and deficient in metals. Here, we present our own 
analysis of the same FEROS spectra, as well as the new 
spectropolarimetric data. As this star varies on such a short 
time scale, we analyze the original intensity spectrum de- 
rived from each single sub-exposure, instead of the co-added 
spectra described in Sect. 2. We limit our current study to 
H i. He I, C II, Si in, and Fe ill. 



line formation cods[E leaving only the assumed circumstellar 
profiles. The subtraction of a single synthetic profile from 
these spectra, assumes that the hydrogen profile does not 
vary from any alternative processes (e.g. pulsation). With 
the new spectra included, we confirm the double S-wave 
variation of circumstellar hydrogen. The equivalent width 
curves for each of the Balmer hydrogen lines show maximum 
emission at phases 0.2 and 0.8, coincident with viewing the 
phases of magnetic extrema. Maximum absorption occurs 
at phases 0.0 and 0.4, coincident with viewing the star's 
magnetic equator. At these phases, we speculate that extra 
absorption may be due to occultation of the star by cir- 
cumstellar material. The secondary photometric minimum 
occurs at the same rotational phase, suggesting a small por- 
tion of the stellar disk is blocked by the co-rotating cloud 
of plasma. The asymmetry of the red and blue shifted emis- 
sion at maximum value confirms that the cloud emission 
measures are not equivalent. 



4.1 Hydrogen 

In Figure 2 of Townsend et al. (2005), the Ha FEROS spec- 
tra are presented in their comparison of RRM model pre- 
dictions with observations. Here, Fig. |4]plots the equivalent 
widths (EW) of the FEROS, ESPaDOnS, and Narval spec- 
tra for Ha, H/3, and H7 as a function of rotational phase. 
Higher Balmer lines show trace amounts of circumstellar 
emission. The dynamical plots of Ha and H/3 (Fig. [5} shows 
the spectra minus a photospheric synthetic spectrum com- 
puted using a 23000 K, log g = 4.0 non-LTE model from 
TLUSTY (Lanz & Hubeny 2007) together with SYNSPEC 



4.2 Helium 

It is well established that a Ori E is a He-strong star with 
patches of helium overabundance surrounded by regions of 
normal or deficient abundance (e.g. Veto 1991; Reiners et 
al. 2000). Fig. [4] shows the equivalent width of the He I 
5876 A and 6678 A lines as a function of rotational phase. 
The dynamical plot in Fig.[S]displays the spectra as a differ- 
ence from the mean profile. Maximum absorption occurs at 

1 http:/ /nova. astro. umd.edu/Synspec43/synspec. html 
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phase 0.4 and 0.9, whereas minimum absorption occurs at 
phases 0.2 and 0.6. The maximum absorption at phase 0.9 
is several times the absorption at phase 0.4. The timing of 
the smaller absorption maximum corresponds with a phase 
of magnetic null, as well as a region where the Ha emis- 
sion is at a relative minimum. The increase in absorption at 
phase 0.4 may then be due to occultation of the star by the 
circumstellar material, as suggested for hydrogen. If this is 
the case, we would expect a similar effect at phase 0.0, how- 
ever the effect should be much smaller than the observed 
absorption maximum at phase 0.9. 

4.3 Metals 

In addition to H and He, we also consider the C II 4267 
A, Si in 4553 A, and Fe ill 5127 A lines. Fig. [J displays 
the equivalent widths as a function of rotational phase for 
these lines. Dynamical plots are shown in Fig. [6] with the 
spectra plotted as the difference from the mean profile. All 
three lines vary similarly, showing an absorption minimum 
at phase 0.9. This feature appears slightly wider in phase 
coverage in silicon, however with the error bars, any differ- 
ence is difficult to detect. Iron shows additional variability, 
although the spectral lines are weak and difficult to normal- 
ize. In the dynamical plot of iron, a feature at phase 0.0 
differs in movement from the rest of the features in these 
plots. It is possible that this excess absorption is due to the 
same circumstellar effect as described for hydrogen and he- 
lium. The minimum absorption of these metal lines at phase 
0.9 coincides with maximum helium absorption. This points 
to decreased metals in helium-strong regions, although it is 
unclear whether there is in fact a second absorption mini- 
mum to correspond to the helium maximum at phase 0.4. 

4.4 Summary 

New spectroscopic data confirm the diverse variability of 
optical spectral lines of a Ori E. Hydrogen emission min- 
ima directly correspond to null longitudinal magnetic field 
measurements. At these rotational phases (0.0 and 0.4), a 
cloud of plasma lies directly in front of the star. The data 
indicate that these clouds are located at the intersection 
between the magnetic and rotational equators. Regions of 
helium enhancement are coincident with regions of metal 
deficiency. The longitudinal field maxima map to phases of 
maximum hydrogen emission. At these phases, we view the 
star's magnetic poles with the plasma clouds in quadrature. 
These new data are in good agreement with not only pre- 
vious observations, but also the overall qualitative plasma 
geometry predicted by the RRM. 



5 INITIAL CONFRONTATION 

As a first exploitation of the current data set, we evalu- 
ate the appropriateness of the magnetic field configuration 
assumed by Townsend et al. (2005). Prior to the present pa- 
per, the longitudinal magnetic field curve for a Ori E could 
be reasonably fit with a sinusoid, indicating a simple dipole 
structure. Townsend et al. (2005) invoked a significant off- 
set of the dipole magnetic field in their application of the 
RRM model to a Ori E in order to explain asymmetries in 
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Figure 7. Longitudinal magnetic field measurements (filled cir- 
cles) for <T Ori E, as well as data reported in Landstreet & Borra 
(1978) and Bohlender et al. (1987) (squares). Each data point is 
plotted with corresponding la error bars. The solid curve is the 
model longitudinal magnetic field curve from the offset dipole con- 
figuration of Townsend et al. (2005). The data are phased using 
the ephemeris of Townsend et al. (2010). 



both the Ha emission spectra and the photometric varia- 
tions. Using the ephemeris of Townsend et al. (2010), the 
offset dipole model fits the historical longitudinal field curve 
with a reduced X 2 of 2.92 (20 de grees of freedom; 22 data 
points - 2 model parameters). The error bars on these data 
are large enough that they are not able to distinguish be- 
tween a large range of magnetic field configurations. When 
the new data are added (Fig. [7]), the model reproduces the 
general character of the observations, but disagrees quanti- 
tatively with the detailed variation. The offset dipole model 
fit to the new measurements computed from LSD profiles 
gives a reduced x 2 value of 92.9 (16 degrees of freedom). 
The new H/3 measurements fit with the same model give a 
reduced % 2 °f H-2 (16 degrees of freedom), while the new 
He line measurements have a reduced \ 2 of 21.6 (34 degrees 
of freedom). 

Since the longitudinal magnetic field is an average quan- 
tity over the stellar disk, an even more stringent test of 
the RRM model magnetic field configuration can be ob- 
tained from the Stokes / and Stokes V profiles, which, by 
Doppler broadening, map variations of physical quantities 
across the surface of the star. Fig.[H]compares the predictions 
of the RRM offset dipole model with observed Stokes I and 
Stokes V profiles of the He I 5876 and 6678 A lines. The syn- 
thetic Stokes parameters were calculated using the forward 
mode of the MDI code InverslO, adopting v sin i=140 km s" 1 
(determined using the current spectra) and i = 75° . In this 
specific model, the He abundance distribution was adopted 
from the Doppler imaging analysis of a Ori E, which will 
be presented in an upcoming paper. Nevertheless, models 
with uniform He distribution applied to the same magnetic 
field parameters produce similar profiles. The magnetic field 
configuration of Townsend et al. (2005) predicts a far more 
extreme phase modulation and much larger amplitudes of 
the Stokes V profiles than are observed. The disagreement 
between the synthetic and observed profiles is especially ap- 
parent for rotational phases 0.318-0.720. 
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Comparison of the RRM model magnetic field configu- 
ration and new spectropolarimetric data in both longitudi- 
nal magnetic field curve and Stokes / and Stokes V profiles 
demonstrate that the specific offset dipole field model in- 
voked by Townsend et al. (2005) is not compatible with the 
new data. The strong offset produces a much larger varia- 
tion of the magnetic field across the visible stellar surface 
than is reflected by the observations, as illustrated by the 
model fit to the Stokes V profiles. These profiles indicate 
that the model requires a very different field topology, fur- 
ther suggested by the periodic variability of the longitudi- 
nal magnetic field. While we cannot definitely rule out the 
possibility of some alternate offset dipole, the asymmetric 
shape of the longitudinal magnetic field curve indicates a 
more complex field, requiring higher-order field components 
to match the shape of the variation. 



6 SUMMARY AND FUTURE WORK 

This paper lays the observational groundwork to re-evaluate 
the RRM model for a Ori E. New spectroscopy is in agree- 
ment with previous observational studies of a Ori E (e.g. 
Pedersen & Thomsen f977; Landstreet & Borra f978; Rein- 
ers et al. 2000). The Ha variability confirms the general 
picture derived from the RRM model of a rapidly rotating 
star in which the magnetic field overpowers the stellar wind, 
allowing plasma to become trapped in a magnetosphere that 
co-rotates with the star. While the RRM model agreed rea- 
sonably well with the previously observed longitudinal mag- 
netic field curve for a Ori E, new, more precise measure- 
ments reveal a substantial variance between the shapes of 
the observed and modeled curves. These results challenge 
the offset dipole form assumed in the Townsend et al. (2005) 
application of the RRM model to a Ori E, and indicate that 
future models of its magnetic field should also include com- 
plex, higher-order components. 

As observations improve, more comprehensive modeling 
is required to explain the intricate physical phenomena of a 
Ori E. We are currently in the process of using MD1 to de- 
rive both abundance maps for strong variable spectral lines 
and magnetic field maps based on several He I lines. The 
circumstellar features observed in line profiles should not 
affect the resulting surface abundance maps, as these fea- 
tures have a much higher radial acceleration than the stellar 
surface (see e.g. Donati & Collier Cameron 1997; Donati et 
al. 1999). These new abundance maps will provide detailed 
information of the surface distribution of several elements, as 
well as provide the information to compute a total synthetic 
photospheric light curve for a Ori E. We plan to use de- 
tailed model atmospheres to calculate flux distributions and 
simulate a light curve, as described by Krticka et al. (2007). 
The magnetic field mapping will offer greater clarity into the 
structure of the magnetic field, however as our data do not 
include linear polarization (Stokes Q and U spectra) data, 
magnetic inversions require an assumption about global field 
topology to reduce ambiguity. In its current state, the RRM 
model can only compute predictions for a dipole magnetic 
field. We anticipate, with these future models, a clearer un- 
derstanding of the remarkable physical phenomena of a Ori 
E. Moreover, this work will provide the groundwork for anal- 
ysis and modeling of similar massive magnetic stars. 
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Figure 8. Observed Stokes / and V profiles of the He I 5876 and He I 6678 A lines and synthetic spectra corresponding to the offset 
dipole configuration of Townsend et al. (2005). There is a clear disagreement between model and data. The phases noted are computed 
using the ephemeris HJD = 2442778^819 + 1"?190842E. The bars at the bottom of the plots give horizontal (1 A) and vertical (5%) 
scales of the profile plots. 
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